* robustness_tpaper_informal_edom.do	JEP	01/26/15
* Does robustness checks that are in the term paper, for informal workers: RIF10 98q4
* Exclude domestic workers

*************************** Change log ********************************** 

/*
	
*/	

clear 
version 13.1

* Project info
cap project, doinfo
if _rc==198 {
	if c(username)=="J16339" loc master "B:\Col_Minimum_Wage" 
	else if c(os)=="Windows" loc master "C:/Users/jorpp/Dropbox (Brown)/Col_Minimum_Wage"
	else loc master "/home/jperez/AA_Minimum_Wage"
}
else {
	local master "`r(pdir)'"
	local doname "`r(dofile)'" 
	* Project calls
	project, uses("`master'/Data/Source/ocupados_mw.dta")
	* Call programs
	project, relies_on("`master'/Code/Est/sample_choice.do")
	project, relies_on("`master'/Code/Est/build_mw_measures.do")
	project, relies_on("`master'/Code/Est/build_treatment.do")
	project, relies_on("`master'/Code/Est/export_rif.do")
	project, relies_on("`master'/Code/Est/confirmedrun.do")
}
*/


* Call necessary programs. I do not call using project to carry the programs around. I do not put these in the ado/personal folder because they are specific to the project. They are not ado files.
do "`master'/Code/Est/sample_choice.do"
do "`master'/Code/Est/build_mw_measures.do"
do "`master'/Code/Est/build_treatment.do"
do "`master'/Code/Est/export_rif.do"
do "`master'/Code/Est/confirmedrun.do"


use "`master'/Data/Source/ocupados_mw.dta", clear

* Choose quantiles
glo qtiles "5 10 15 20 25 30 50 75 90"
* Number of bootstraps
glo boot=10

* RIF 4 - Tables numbered 4 in term paper

* Choose sample
mysample rif10sample , hmin(30) hmax(50) sex(0) edadm(65) wage(salario) informal
tab rif10sample
tab tipo_trabajador if rif10sample
tab year if rif10sample
tab horas_semana if rif10sample 
tab sex if rif10sample
tab afiliado_salud if rif10sample
* Build mw measures
mw_measures rif10sample , stub(rif10)
* Build treatment
mytreatment, posts("tq(1998q4)") mw(bite_rif10) wage(salario_hora_real)

* RIF regressions
* Create estimation results folder if it doesn't already exist
cap mkdir "`master'/Ster/rif10_new"

* I found that these do not replicate what was initially in the paper exactly, because there were some observations that were being mistakenly included. Coefs are pretty close anyway.

* This is to be able to build the project a second time without re-running the regressions if the .ster files are already there


foreach x in $qtiles {
	di " Percentile `x'"
	loc fw: word 1 of $qtiles
	loc count: word count $qtiles
	di `count'
	loc lw: word `count' of $qtiles
	if `x'==`lw' loc reg="register"
	else loc reg ""
	if `x'==`fw' loc r="replace"
	else loc r="append"
	* Panel A. Employment
	*confirmedrun, `reg' controls("") name("`master'/Ster/rif10_new/reg_emp_tq(1998q4)_bite_salario_hora_real_50_edom") rif(rif10A) p(`x') r(`r') sample(rif10sample)
	* cgmrif2 lw treat i.ccode i.time emp if rif10sample , hyp(0) boot($boot) cluster(ccode) x(treat) i(ccode) t(time) q(`x') rif(rif10A`x')
	* exportrif, name(`master'/Ster/rif10_new/reg_emp_tq(1998q4)_bite_salario_hora_real_50_edom) p(`x') r(`r')
	* Panel B. Employment and Bartik Price Variable
	confirmedrun, `reg' controls(bartikp1) name("`master'/Ster/rif10_new/reg_emp_bp_tq(1998q4)_bite_salario_hora_real_50_edom") rif(rif10B) p(`x') r(`r') sample(rif10sample)
	* Panel C. Employment and Bartik Quantity Variable
	*confirmedrun, `reg' controls(bartikq) name("`master'/Ster/rif10_new/reg_emp_bq_tq(1998q4)_bite_salario_hora_real_50_edom") rif(rif10C) p(`x') r(`r') sample(rif10sample)
	* Panel D. Employment and City Trends
	*confirmedrun, `reg' controls(i.ccode#c.trend) name("`master'/Ster/rif10_new/reg_sst_emp_tq(1998q4)_bite_salario_hora_real_50_edom") rif(rif10D) p(`x') r(`r') sample(rif10sample)
	* Panel E. Employment, Bartik Price Variable and City Trends
	*confirmedrun, `reg' controls(i.ccode#c.trend bartikp1) name("`master'/Ster/rif10_new/reg_sst_emp_bp_tq(1998q4)_bite_salario_hora_real_50_edom") rif(rif10E) p(`x') r(`r') sample(rif10sample)
	* Panel F. Employment, Bartik Quantity Variable and City Trends
	*confirmedrun, `reg' controls(i.ccode#c.trend bartikq) name("`master'/Ster/rif10_new/reg_sst_emp_bq_tq(1998q4)_bite_salario_hora_real_50_edom") rif(rif10F) p(`x') r(`r') sample(rif10sample)
}
